Задача о погоне
Чемоданова Ангелина Александровна
Российский университет дружбы народов имени Патриса Лумумбы, Москва, Россия
26 февраля 2025
Построить математическую модель для выбора правильной стратегии при решении примера задаче о погоне.
На море в тумане катер береговой охраны преследует лодку браконьеров. Через определенный промежуток времени туман рассеивается, и лодка обнаруживается на расстоянии 25 км от катера. Затем лодка снова скрывается в тумане и уходит прямолинейно в неизвестном направлении. Известно, что скорость катера в 5,1 раза больше скорости браконьерской лодки.
Записать уравнение, описывающее движение катера, с начальными условиями для двух случаев (в зависимости от расположения катера относительно лодки в начальный момент времени).
Построить траекторию движения катера и лодки для двух случаев.
Найти точку пересечения траектории катера и лодки
Чтобы найти расстояние x (расстояние после которого катер начнет двигаться вокруг полюса), необходимо составить простое уравнение. Пусть через время t катер и лодка окажутся на одном расстоянии от полюса. За это время лодка пройдет x , а катер k − x (или k + x, в зависимости от начального положения катера относительно полюса). Время, за которое они пройдут это расстояние, вычисляется как $\dfrac{x}{v}$ или $\dfrac{k-x}{5.1v}$ (во втором случае $\dfrac{k+x}{5.1v}$). Так как время одно и то же, то эти величины одинаковы. Тогда неизвестное расстояние можно найти из следующего уравнения: $$ \dfrac{x}{v} = \dfrac{k-x}{5.1v} \text{ -- в первом случае} $$ $$ \dfrac{x}{v} = \dfrac{k+x}{5.1v} \text{ -- во втором} $$ Отсюда мы найдем два значения $x_1 = \dfrac{25}{6,1}$ и $x_2 = \dfrac{25}{4,1}$, задачу будем решать для двух случаев.
После того, как катер береговой охраны окажется на одном расстоянии от полюса, что и лодка, он должен сменить прямолинейную траекторию и начать двигаться вокруг полюса, удаляясь от него со скоростью лодки v. Для этого скорость катера раскладываем на две составляющие: vr - радиальная скорость и - vτ тангенциальная скорость. Радиальная скорость - это скорость, с которой катер удаляется от полюса, $v_r = \dfrac{dr}{dt}$. Нам нужно, чтобы эта скорость была равна скорости лодки, поэтому полагаем $\dfrac{dr}{dt} = v$.
Тангенциальная скорость – это линейная скорость вращения катера относительно полюса. Она равна произведению угловой скорости $\dfrac{d \theta}{dt}$ на радиус r, $r \dfrac{d \theta}{dt}$.
Получаем:
$$v_{\tau} = \sqrt{26.01v^2-v^2} = \sqrt{25.01}v$$
Из чего можно вывести:
$$ r\dfrac{d \theta}{dt} = \sqrt{25.01}v $$
Решение исходной задачи сводится к решению системы из двух дифференциальных уравнений:
$$\begin{cases} &\dfrac{dr}{dt} = v\\ &r\dfrac{d \theta}{dt} = \sqrt{25.01}v \end{cases}$$
С начальными условиями для первого случая:
$$\begin{cases} &{\theta}_0 = 0\\ \tag{1} &r_0 = \dfrac{25}{6.1} \end{cases}$$
Или для второго:
$$\begin{cases} &{\theta}_0 = -\pi\\ \tag{2} &r_0 = \dfrac{25}{4.1} \end{cases}$$
Исключая из полученной системы производную по t, можно перейти к следующему уравнению: $$ \dfrac{dr}{d \theta} = \dfrac{r}{\sqrt{25.01}} $$
Начальные условия остаются прежними. Решив это уравнение, мы получим траекторию движения катера в полярных координатах.
# функция, описывающая движение лодки браконьеров
x(t) = tan(fi)*t;
# функция, описывающая движение катера береговой охраны
f(r, p, t) = r/sqrt(25.01)
# постановка проблемы и решение ДУ для 1 случая
prob = ODEProblem(f, r0, theta0)
sol = solve(prob, saveat = 0.01)
# отрисовка траектории движения катера
plot(sol.t, sol.u, proj=:polar, lims=(0, 15), label = "Траектория движения катера")Получим график движения катера в 1 случае:
Добавим в код:
Получим график движения катера и траекторию движения лодки в 1 случае:
Теперь рассчитаем точку пересечения траектории катера и лодки.
# точное решение ДУ, описывающего движение катера береговой охраны
y(x)=(250*exp((10*x)/(sqrt(2501))))/(61)
# подставим в точное решение угол, под которым движется лодка браконьеров для нахождения точки пересечения
y(fi + pi)Получим следующее значение: 12.304002757914663.
При оценки точки на глаз можно получить похожее значение:
Перейдем к решению второго случая.
Получим график движения катера во 2 случае:
Добавим в код:
Получим график движения катера и траекторию движения лодки во 2 случае:
Теперь рассчитаем точку пересечения траектории катера и лодки.
Получим следующее значение: 9.767236102657977. При оценки точки на глаз можно получить похожее значение:
Мы построили математическую модель для выбора правильной стратегии при решении примера задаче о погоне.